Near-field heat transfer between a nanoparticle and a rough surface 
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In this work we focus on the surface roughness correction to the near-field radiative heat transfer 
between a nanoparticle and a material with a rough surface utilizing a direct perturbation theory 
up to second order in the surface profile. We discuss the different distance regimes for the local 
density of states above the rough material and the heat flux analytically and numerically. We show 
that the heat transfer rate is larger than that corresponding to a flat surface at short distances. At 
larger distances it can become smaller due to surface polariton scattering by the rough surface. For 
distances much smaller than the correlation length of the surface profile, we show that the results 
converge to a proximity approximation, whereas in the opposite limit the rough surface can be 
replaced by an equivalent surface layer. 



I. INTRODUCTION 

Recently, the near-field radiative heat transfer has at- 
tracted a lot of theoretical and experimental attention^— . 
It was predicted theoretically 4 - and shown experimen- 
tally^— that the heat flux for distances much smaller 
than the thermal wavelength Ath = hc/{k-QT) can be 
much greater than that predicted by Planck's law, where 
h is Planck's constant, k& is Boltzmann's constant, c is 
the velocity of light in vacuum and T the temperature. 
This unusual property might for example be exploited for 
thermo-photovoltaics^— and near-field scanning thermal 
microscop y 13 ' 14 . 

It is common knowledge that the radiative properties 
of a material depend not only on the material parameters 
but also on the surface roughness^ 5 .. While the effect on 
far-field properties has been widely studied^, the impact 
of surface roughness on near-field heat transfer has not 
been considered so far. From the experimentalist's point 
of view, at least an estimate of the surface roughness cor- 
rection is desirable, since one is confronted with a certain 
degree of surface roughness in all experimental set-ups. 

In this work we will study the near-field heat trans- 
fer between a nanoparticle considered to be an electric 
dipole and a material with a rough surface. Within this 
model the change of the local density of states (LDOS) 
above the material due to the surface roughness com- 
pletely causes the change in the heat flux. Since, the 
near-field heat transfer between two semi-infinite bodies 
is also largely determined by the LDOS, we think that 
our results are not only restricted to the here discussed 
geometry, but can also be utilized to get a rough estimate 
of the impact of surface roughness for configurations used 
in recent experimental setups 5-8 . In addition, the LDOS 
also plays a fundamental role for other physical phenom- 
ena as it determines, for instance, the lifetime of atoms 
and molecules near a surface, so that the here given re- 
sults for the LDOS have a wider range of application. 

The paper is organized as follows: In Sec. II we give 
a short description of the dipole model of near-field heat 
transfer. In Sec. Ill we introduce the perturbation result 
for the mean LDOS, the key quantity for understanding 



the roughness correction to the mean heat flux which is 
itself discussed in Sec. IV. In Sec. V we derive approxima- 
tions for the small and the large distance regime. Finally, 
in Sec. VI we study the roughness correction to the LDOS 
numerically 



II. RADIATIVE HEAT TRANSFER 

We consider the situation depicted in Fig. [TJ A 
nanoparticle with a polarizability a in local thermal equi- 
librium at temperature Tp is placed at rp near a dielectric 
half-space with a given permittivity e. As discussed in 
Refsi 17 ' 18 the multiple scattering between the nanoparti- 
cle and the surface can be neglected. We assume that this 
dielectric is in local thermal equilibrium at a temperature 
Tb 7^ Tp. The interface separating the dielectric from the 
vacuum is described by the surface profile S(x) — z = 
with x = (x, y). 
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Figure 1: Sketch of the configuration considered here. 

Now, as far as the radius R of the nanoparticle is 
smaller than the thermal wavelength Ath and the dis- 
tance d between the dielectric body and the particle can 
be assumed to be large, i.e., A t h ^> R and d ^> R, the 
energy transfer rate between the particle and the dielec- 
tric body due to radiation can be expressed within the 
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dipole model as (see Refsj 2 -^ and references therein) 
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where /3 = (fceT) -1 , &b is Boltzmann's constant and 
{2irh) is Planck's constant. Here, the spectral power 
absorbed by the nanoparticl o 2 ' 3 ' 19 is given by the term 
a"(uj)D E (uj, rp)0(w, Tb), i.e., it is proportional to the 
imaginary part of the polarizability a" of the particle and 
proportional to the electric energy density which is given 
by the product D E (io, rp)0(w, Tb) of the electric local 
density of states (LDOS) above the dielectric medium 
and the mean energy of a harmonic oscillator Q(ui,T). 
On the other hand the power emitted by the particle 
and absorbed within the bulk medium is proportional to 
a"(uj)D E (uj, rp)0(w, Tp). We point out that the above 
given formula has to be augmented by its magnetic coun- 
terpart when considering metallic materials as discussed 
in Ref3. 

Within this work we will use the expression of the 
LDOS 

D E ( Wl r) = ^ImTrG E (r,r) (3) 

7TCT 

as defined in Ref—. Using Eq. ([3]) in Eq. (JlJ corresponds 
to a situation where the bulk and its surrounding are 
at ambient temperature Tb, whereas the nanoparticle is 
heated or cooled with respect to Tb. 



III. LOCAL DENSITY OF STATES 

A. Stochastic surface roughness 

In this work, we concentrate on the special case of 
a stochastic surface profile S described as a stochastic 
Gaussian process with mean value and correlation func- 
tion given by 



<S(x)) p = 0, 
(S(x)S(x')) p =8 2 W(\x-x'\). 



(4) 
(5) 



The brackets (o) p stands for the average over an ensemble 
of realizations of the surface profile 5(x); 8 is the rms 
height of the surface profile. The correlation function 
W(|x — x'|) is here assumed to be given by a Gaussian 



W(\x-x'\) = e~ 



(6) 



introducing the transverse correlation length a. 

For the Fourier component S(k) of the surface profile 
function one obtains 

(S(k)> p = (7) 
(S(k)S(k')) p = (2tt) 2 8 2 8{k + k') 9 (k) (8) 



with Dirac's delta function S(k + k') and the surface 
roughness power spectrum 



9 (k) = d 2 xW(\x\)e- 
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B. Perturbation expansion of the LDOS 



In order to determine the perturbation expansion of 
the LDOS in Eq. ([3]), we expand the Green's dyadic G E 
with respect to the surface profile up to second order—. 
We follow Ref— , where one can find a procedure for 
determining the first-order Green's dyadic. The explicit 
second order form is given in appendix A. A detailed 
discussion of the validity of the perturbation theory is 
given in Ref^. In summary, the perturbation ansatz is 
valid as far as the surface height 8 is the smallest length- 
scale of the problem, i.e., 8 <C min{z, a, Ath}- 

Inserting the Green's function up to second order given 
in Eq. ([Al) into Eq. ©, we find for the LDOS up to 
second order after ensemble average: 
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where /i. 



1 and h r 



(2 k 2 - fc 2 )/fc 2 . The first term 
yields the propagating mode contribution, i.e., for k < fco 
with 7 r = yj /lq — k 2 purely real, whereas the second term 
gives the contribution due to evanescent waves, i.e., k > 
ko with 7 r = i'y purely imaginary. The mean reflection 



coefficient up to second order 
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can be written 



as a sum of the usual Fresnel coefficient r s / p and a surface 
roughness correction 



(^/p" (2) >P = Vp-2i7 r P°/p) 2 Af s / P , (11) 



where 
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and 
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with 7t = y/kge — k 2 . The so-called proper self-energy 
M s / p is defined in Appendix |B] Furthermore, one can 
write the expression for the proper self energy M s / p as a 
sum of two contributions M s / p o originating from terms 
due to a second-order scattering process and M s / p l orig- 
inating from terms due to two successive first-order scat- 
tering processes (see Fig. [2]). 



IV. RADIATIVE HEAT TRANSFER BETWEEN 
A NANOPARTICLE AND A ROUGH SURFACE 

With the above given relations we can study the radia- 
tive heat transfer between a spherical nanoparticle and 
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Figure 2: Simple sketch of the scattering processes due to 
Ms/p,o ( a ) an d due to M s / p> i (b). After averaging the trans- 
lational symmetry is restored so that k is the same before and 
after the scattering with the rough surface. Mathematically, 
this property follows directly from the statistical properties of 
the gaussian surface roughness, i.e., (S' 2 ' (#4— k'))p k S(k-k') 
and (S w (k - k')S (1) (k' - k")>p oc 6(k - k"). 



a semi-infinite dielectric body with a rough surface for- 
mally given by Eq. (flj setting the nanoparticle at posi- 
tion rp = (0, 0, d). For describing the absorptivity of the 
nanoparticle with radius R we utilize the polarizability 
given as 



a(u>) = AirR 



3 e ( w ) 



e(«) 



(13) 



Here, we employ the material properties for SiC from 
RefiS4 for numerical evaluation of Eq. ([T]) using Eq. (|10p 
and determine the surface roughness correction defined 
by 



/p(0)-(2)\ _p(0) 



(14) 



We have checked that P (0) gives the same result as in 
Refii when choosing the same radius. (Note that AP 
itself does not depend on R, since (P(°' - ' 2 ') p oc R 3 and 
P(°) oc R 3 .) 




10' 10" 
d (meters) 

Figure 3: (Color online) Plot of the modulus of AP for evanes- 
cent modes over the distance for SiC setting Tb = OK and 
Tp = 300 K for a rough surface with S — 5 nm and a = 200 nm. 
The red part of the curve indicates positive and the blue one 
negative values. The PA and the LDA are included for com- 
parison. 



Let us first turn to the distance dependence of the 
heat flux. A plot of AP over the distance considering 
only evanescent modes, with 8 — 5 nm and a = 200 nm 
is shown in Fig. [3J The temperature of the dielectric 
is assumed to be OK, whereas the temperature of the 
nanoparticle is set to 300 K. From formula (PQ) it is clear 
that (apart from a sign for P) one gets the same result 
for AP when interchanging the temperatures. As will be 
shown later, the proximity approximation (PA) 



(P)l A = P<°) 1 + 6 



.S 2 
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(15) 



can be derived in the small distance regime with ()Cfi< 
a. Note that AP is always positive and independent from 
T in that limit even if e depends on T. In the opposite 
limit with d a, we can also derive a large distance 
approximation (LDA) as will be shown later (see Eq. 
and 
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Figure 4: (Color online) As Fig. [3] but for evanescent and 
propagating modes. 

It can be seen, that AP converges to the approxima- 
tions for large and small distances. For distances slightly 
greater than the small distance regime well described by 
the PA result the surface roughness correction becomes 
much greater than predicted by the PA. For distances 
between 1—10 /im it can be seen that AP becomes neg- 
ative. Now, in this distance regime already the prop- 
agating modes start to dominate the heat transfer. It 
turns out that for the propagating modes the surface 
roughness correction is very small compared to that of 
the evanescent modes. It follows that AP for evanescent 
and propagating modes becomes also very small in the 
distance regime 1 — 10 fxm as is illustrated in Fig.HJ It can 
also be seen that due to the competition of the roughness 
correction of the evanescent and propagating modes, the 
overall correction to the heat transfer AP becomes more 
oscillatory in that distance regime. Nevertheless, correc- 
tions on the order of ten percent can arise for distances 
smaller than 100 nm, but one has to keep in mind that 
the theory used here is only applicable for d 3> S and 
d 3> R- Therefore, using 5 = 5 nm and R = 5 nm it can 
be estimated that the theory is only valid for distances d 
greater than ss 50 nm. 
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Figure 5: (Color online) As Fig. [3] using a — 100 nm, 200 nm 
and 500 nm. 

Since wc have determined P pcrturbatively with re- 
spect to the surface profile it is clear that to second order 
AP is proportional to S 2 . On the other hand, the depen- 
dence of AP on the correlation lenght a is not obvious. 
Wc find that AP is approximately proportional to a -1 in 
the large distance regime (for the here used parameters) 
what might be seen in Fig. [5] considering only evanescent 
modes for a = 100 nm, 200 nm and 500 nm. 

In order to understand the roughness correction to the 
heat flux, a deeper understanding of the roughness cor- 
rection to the LDOS is required. Therefore, we will dis- 
cuss the LDOS in more detail in the following. 

V. APPROXIMATIONS OF THE LDOS FOR 
SMALL AND LARGE DISTANCES 

The key quantities for the LDOS are the reflection co- 
efficients in Eq. (jllj) . Here we will derive some approxi- 
mations of the proper self energy M s / p first. From these 
expressions one can get the corresponding approximation 
for the LDOS by using Eq. (JO) and Eq. (fTI). 

Before we derive the small and large distance approx- 
imation for the proper self energy M s / p we first imple- 
ment the following approximation: For frequencies rel- 
evant at room temperature, i.e., u> f=s 10 14 s _1 , we have 
\ko^/ea/2\ <C 1, i.e., a is much smaller than the skindepth 
d s , as far as the correlation length a is small enough. 
Considering SiC, this relation is well fullfillcd for most 
frequencies within the range 3.7 ■ 10 13 — 2 • 10 14 s _1 for 
values of a smaller than 500 nm. Therefore, we concen- 
trate here on the limit a <C d s only (see appendix [C]) , 
while one can determine the opposite limit a>d s with 
a similar procedure. For the latter limit we just state the 
results in Appendix [Dl 

A. small distance regime (i5 < z < a) 

In the limit of small distances z it is seen from Eq. (fTU)) 
that the main contribution comes from wavevectors n rj 
1/z. In the regime a ^> z (na S> 1), the terms M s / Pjl 
become negligible compared to M B / p . This is due to 



the fact that for na 1 the proper self-energy con- 
tributions |M s / p l | are proportional to 1/k or k, resp., 
because the Bessel functions are in this limit approxi- 
mated by exp(£ 2 /8)/£. On the other hand |M s / p | is 
for large wave vectors proportional to k or k 3 , so that 
|-^s/p,i|/|-^s/p,ol 1/k 2 . This means, that in the near- 
field regime the surface roughness correction solely stems 
from A/g/p.o m Eq. (|B10[) and (|BL2|) . i.e., from that term 
which originates from one second-order scattering pro- 
cess (see also Fig. 0). Since the term M s / Pil originating 
from terms due to two successive first-order scattering 
processes including processes involving the excitation of 
surface modes with (k, u>) followed by scattering into an- 
other surface mode («', ui) as an intermediate state which 
is then scattered back into the initial state with (k, u>) 
becomes negligible, we can formally conclude that these 
processes are irrelevant for z <C a. This is very intuitive, 
since the wavelength of the surface modes becomes much 
smaller than the correlation length. Within this limit the 
excited surface modes can follow the perturbed surface 
adiabatically without being scattered. Hence, for z <C a 
it suffices to consider M s / P:0 only. Inserting M s / Pi o into 
the mean reflection coefficient in Eq. (fTTj) gives for the 
evanescent near-field regime with k»I;o 

C (2) )p^s/p(H2M 2 ), (16) 

for > 1 or z < a, utilizing the quasistatic approxi- 
mations for the reflection coefficients 

r s = ^(e-l) and r p = ^j. (17) 

From the above relations we can infer that the surface 
roughness correction for very small distances does not 
depend on a. The expressions in Eq. (|16[) can also be 
derived for a d s . 

By inserting Eqs. (fT()]) and (jTTJ) into the corresponding 
formula for the LDOS in the quasistatic limit to second- 
order perturbation theory 

(18) 

which follows from Eq. fj 10[) when assuming k > Ico, we 
find 

(P E (z))P A «P E(0) (l + 6^). (19) 

From this result the approximation for the radiative heat 
transfer in Eq. (| L5[) follows easily. 

Now, exactly the same result can be obtained with the 
so-called proximity approximation (PA), which holds in 
the quasistatic limit for z <ti a and z 3> 5 which was for 
example used to determine the surface roughness contri- 
bution to the Casimir forc o 26 i 27 and has recently been 
employed to the near-field radiative heat transfer—. It 
amounts to replace the rough surface by a horizontal 
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plane at a random height z = 5(x) followed by the en- 
semble average. Hence, the LDOS reads: 
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D E{0) S 2 + O(4). 



(20) 



Utilizing the expression for the LDOS in the quasistatic 
limitSi 
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one retrieves the PA in Eq. (|19p when considering terms 
up to second order only. Since the zeroth-order LDOS 
is proportional to l/z 3 , it is clear that the contributions 
for planes at distances smaller than z will give a larger 
value than those at distances larger than z so that af- 
ter averaging the overall roughness correction is positive. 
Regarding only second-order terms it is also clear that 
this correction is proportional to 5 2 . 



Therefore, for frequencies near the surface resonance fre- 
quency and for distances z 3> a the proper self energy 
M s / p in the expressions for the reflection coefficient in 
Eq. (JTTJ) can be replaced by Eqs. an d ((Ml) yielding 
the large distance approximation (LDA). 

As far as one considers correlation lengths a much 
smaller than the wavelength inside the medium X/\y/e\, 
one can apply the concept of homogenization 3 -^ by replac- 
ing the surface roughness by an equivalent surface layer 
with an effective permittivity. In the near-field regime, 
this condition is fullfilled if na <C 1 or a <C z, resp. Fol- 
lowing the approach of Rahman and Maradudi n 32 ' 33 we 
replace the rough surface in the limit of large distances 
with z 3> a and for a <C d s by a thin equivalent sur- 
face layer as depicted in Fig. [6] ranging from z = —ah 
to (1 — a)L with a € [0, 1]. The permittivity £l of that 
layer is considered to be 0.5(e+ 1), i.e., the average of the 
permittivity of vacuum and the dielectric medium. The 
reflection coefficient for such a layered geometry with 
£i = e, £2 = £Lj and £3 = 1 is in the quasistatic limit 
K > fco 



B. 



large distance regime (A t h >z>a and a <C d s ) 



When considering the evanescent contribution to the 
heat transfer in the large distance regime with d>a one 
is facing a situation as sketched in Fig. [6] (a). At room 
temperature Ath ~ 10 /an, so that the evanescent regime 
is restricted to d < Ath = 10 /xm. Therefore, this large 
distance regime is only applicable for surface profiles with 
a correlation length a much smaller than 10 /im and a rms 
i « s. The same is true for the LDOS at z > a. On 
the other hand, for very small temperatures of about 5 K, 
the thermal wavelength is about A t h ~ 0.46 mm. Hence, 
for low temperature experiments as the measurement of 
spin-flip lifetimes in Rcfi 29 ' 30 , this large distance regime 
can be applied to a much wider range of surface roughness 
parameters, distances and nanoparticle radii. 

Considering the large distance limit na <C 1 in 
Eq. (|C1[) . we find for the lowest nonvanishing order 
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and with Eq. JC4 
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(23) 



Obviously, M s / p l are proportional to a -1 . Therefore, for 
large distances for which |M s / P) jJ ^> |M s / p0 | is fullfilled, 
the second-order correction to the LDOS will be inversely 
proportional to the correlation length if a <C d s . The con- 
dition |M s / p l | 3> |M s / p | is fullfilled for the p-polarized 
modes in the here given limit na -C 1, whereas for the s- 
polarized modes, this is only true if 27 t a(e+ 1) <C \/tt. In 
addition, the approximate expressions in Eqs. (|22j) and 
(|23p have the denominator e + 1 indicating a strong con- 
tribution for surface resonances with e' = — 1 and e" « 0. 
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Figure 6: Sketch of the physical situation encountered in the 
large distance regime for the heat transfer (a) and the replace- 
ment of the rough surface by an equivalent layer for surface 
roughness (b) to determine the LDOS in that regime. 

Considering first p-polarised modes we insert the elec- 
trostatic expressions for the Fresnel coefficients 



„32 



(-L + 1 



and 
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(25) 



and expand the reflection coefficient for the layered sys- 
tem with respect to the thickness L, which is thought to 
be small, yielding 



1 



kL 



( e + l) 3 



[~a(l + e) 2 + e 2 + e+l]. (26) 



Now, we want to compare Ar 



r p with the 



result of the LDA. Inserting Eq. (|23)l into the expression 
for the mean reflection coefficients in Eq. (|11[) gives for 
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the quasistatic limit (k 3> ko) 
A(rW-( 2 )) p 



„(0)-(2) 



v^F n5 2 (e - l) 2 
T a (e + 1) 3 



(2,-1). 
(27) 

In order to get a relation between L and a, we compare 
the leading term in e for 8r^ of the layered geometry with 
the surface roughness correction. This gives 



(l-a)L = 



V^ 2 



(28) 



Obviously, the choice of the parameters a and L is 
ambiguous. This ambiguity is resolved in Refi^ 2 - when 
further considering the transmitted field component for 
the equivalent layer and the corresponding perturba- 
tive result. From this, Rahman and Maradudin find 
L = 3y/ir5 2 /a. Using this result in Eq. gives a = 2/3 
as was also found by Rahman and Maradudin in RcfiS 
Hence, the surface roughness can be mimicked by a small 
layer with a thickness of order S 2 /a which is slightly 
shifted into the vacuum region with respect to z = 0, 
so that the effective distance from the surface becomes 
z — L/2>. It follows, that for distances z 3> a the evanes- 
cent part of the LDOS and of the heat transfer, which 
have a monotonous decay with z, will be slightly bigger 
than that of a flat surface and proportional to a~ 1 . 

The same considerations can be made for the s- 
polarized modes yielding the results 
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in the quasistatic regime, which can be compared to the 
corresponding correction to the mean reflection coeffi- 
cient (by inserting Eq. ([52} into Eq. (fTTjll 



A(^ 2) ) P «-(M) 2 ^-i^-^. 

2 na 4 (e + 1) 

Comparing the leading order term of Sr^ and S( 
in e gives 
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Again, the choice of a and L is ambiguous. Using again 
the result for L from Refi^ 2 - gives a = 29/12. Obviously, 
a > 1 so that it is not possible to describe the effect of 
surface roughness for s-polarized modes within this model 
when choosing the same thickness L as for p-polarized 
modes. In order to get the correct L, it is necessary to 
consider the transmitted field contribution. Therefore, 
here it can only be concluded that for the s-polarized 
part we can also use the equivalent layer model with L 
of the order S 2 /a. 

As it might be shown in the next section, in the large 
distance regime, the roughness plays a significant role 
only at the surface polariton resonance, i.e., for u) s such 
that Re(e(w s )) + 1 w and Im(e(w s )) w 0. This is 



clear when noting the resonant denominator l/(e+ 1) 3 in 
Eq. (|2~T|) . Hence, we expect a significant modification of 
the LDOS and therefore also of the lifetime of an emitter 
whose frequency is close to lj s . On the other hand, for 
the near-field heat transfer one has to sum up all contri- 
butions over the spectrum close to the thermal frequency 
so that one can in general expect that the correction is 
in this case small. 

Finally, we note that, although the physics is differ- 
ent, the effect of roughness on heat transfer is the same 
for evanescent and propagating waves. The roughness 
can be modeled by an effective layer with intermediate 
optical properties. In both cases, it results in a larger 
transmission. 



VI. NUMERICAL RESULTS FOR THE LDOS 
A. Propagating modes 

First, we discuss some numerical results for the propa- 
gating modes using the material parameters for SiC from 
Refi 2 ^. In Fig. [3 we show a plot of the deviation of re- 
flectivity from Eq. (fTTj) defined as 



Ai? s/p = 100 
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for s- and p-polarization for k < fco choosing a rms of 5 = 
5 nm and a correlation length of a = 200 nm. Further- 
more, the plot is restricted to frequencies around the sur- 
face phononon polariton resonance u> s = 1.787 • 10 14 s _1 
within the reststrahlcnband, i.e., oj t < u> < u>i, where 
U)l = 1.827 • lO^s" 1 and cj t = 1.495 • 10 14 s _1 are the 
frequencies of the LO and TO phonons in SiC, resp. It 
can be seen that the correction is small and negative in 
this frequency range. For frequencies around the sur- 
face phonon frequency one finds a relatively large neg- 
ative correction (but still smaller than one percent) for 
both polarisations. Hence, the roughness correction to 
the LDOS for propagating modes will be very small as 
well. 

The surface roughness allows for coupling of incident 
propagating waves with surface polaritons 3 ^— . Hence, 
for (w, k) for which the conditions for coupling of surface 
polaritons with propagating modes are met, the reflectiv- 
ity will decrease, i.e., Ai? p is negative, since only a small 
fraction of the excited surface mode will be reradiated. 
This surface wave mediated decrease of the reflectivity 
can be seen for both s- and p-polarization. 

The possibility of exciting surface polaritons with s- 
polarizcd light needs some explanation, since after tak- 
ing the ensemble average, the rough surface has the same 
symmetries as a flat surface for which surface polari- 
tons can be excited with p-polarised light only. On the 
other hand, for surfaces with a grating which breaks 
the translational and rotational symmetry, it is known 
tha t 37 ' 38 s-polarized waves can excite surface polaritons. 



7 



(a) 




Figure 7: (Color online) Plot of Ai? s/p as defined in Eq. ((52)) 
for (a) s- and (b) p-polarized modes using 5 = 5 nm and a = 
200 nm. 




Figure 8: (Color online) Illustration of the excitation of sur- 
face polaritons by s-polarized light in the K-plane (dashed line 
= dispersion relation of the surface phonon polariton, gray cir- 
cle = propagating waves with k < ko): An s-polarized wave 
with Ki = (hi, 0) and an electric field in y direction. Due to a 
first scattering process labeled as 1 with the rough surface the 
surface power spectrum provides the necessary extra momen- 
tum (k[ — K sp ) to match the phase with the surface polariton 
with wave vector K sp . The incoming s-polarized wave has an 
electric field component in the direction of n Bp so that it can 
excite the surface polariton. Due to the second scattering 
process 2 the surface power spectrum again provides the nec- 



essary extra momentum — (m 
s-polarized wave with k = n;. 



K sp ) resulting in a scattered 



In this case, the grating provides the extra momentum for 
matching the phase of the incoming light with that of the 
surface polariton. Additionally, the incoming wave must 
have an electric field vector component parallel to the 
surface phonon polariton wave vector. When these two 
conditions are me t 37 ' 38 , an s-polarized wave can excite 



surface polaritons. Now, for a rough surface the transla- 
tional and rotational symmetry are also broken so that 
an s-polarized wave should be able to excite surface po- 
laritons. An example of such a scattering process is il- 
lustrated and discussed in Fig. [5] 



B. Evanescent modes 

Now, we turn to the evanescent modes. According to 
Eq. (fTU|) . we are interested in Im(r). In order to discuss 
the change in Im(r) we define 



AIm(r s/p ) = 100 



V\ s/p 



Im(r s / F 



(33) 



In Fig. [3] we show a plot of this quantity for s- and p- 
polarized modes using again the parameters (5 = 5 nm 
and a = 200 nm. We show here only the results for na < 
4, i.e., in that region where the mean reflection coefficient 
is dominated by M s / p l . For s-polarizcd modes one can 
see that Im(r s ) is increased up to 60 percent for values of 
(lo, k) coinciding with the dispersion relation of surface 
phonons. The underlying mechanism is the same as for 
the propagating modes depicted in Fig. [5J but with the 
difference that Kj is greater than kg. 



(a) 




Figure 9: (Color online) Plot of Im(r S/ / p ) as defined in Eq. (|33[) 
for (a) s- and (b) p-polarized modes using 8 = 5 nm and a = 
200 nm. 

On the other hand, for p-polarized modes Im(r p ) is 
decreased by about 20 percent for values of (to, k) coin- 
ciding with the dispersion relation of the surface phonons. 
Around na w 1 one also finds a large increase of about 



8 



15 percent for frequencies slighty below and above the 
surface phonon frequency. This can be easily under- 
stood by looking at Fig. ITUl showing a plot of Im(r p ) and 

Im^rp ^ ^)) for na = 1 and for frequencies near the sur- 
face phonon frequency. As can be observed, the disper- 
sion relation is broadened due to roughness induced scat- 
tering of the surface phonons into other surface phonon 
states. For a slightly rougher surface with S = 10 nm 
the broadening becomes a splitting of the surface polari- 
ton dispersio n 33 ' 40 , which can be easily understood from 
the fact that the rough surface acts as a thin layer as 
discussed above. It follows, that the quantity AIm(r s / p ) 
has negative values for frequencies around the surface 
phonon frequency and positive values slightly below and 
above that surface phonon frequency as can be seen in 
Fig. 13 We remark that for 5 = 10 nm the direct pertur- 
bation theory, while qualitatively correct, starts to give 
quantitatively wrong results for frequencies near the sur- 
face resonance frequency 3 ^. 




o»/10"s-' 



Figure 10: Plot of Im(r p ) and Im((r p ' } p ) for na — 1 using 
5 = 5 nm and a = 200 nm. 



C. Spectrum of the LDOS 

From the above discussion it is clear that the roughness 
correction of the LDOS defined as 



AD E = 100 



{D E(0)-(2) )p _ D E(0) 



D E 



(0) 



(34) 



will have positive and negative contributions for different 
frequencies and distances. To make this point clear, we 
-)E a „j n E _ / d e(°)-( 2 ) 



show some plots of AD" and D k 



for 



frequencies ranging from 10 13 s 1 to 2.5 ■ 10 14 s 1 . 

In Fig. [TT] we plotted D E and AD E for propagating 
and evanescent modes and for propagating modes only 
for a distance of z — 5 fim. It is seen that D E is domi- 
nated by the propagating modes and shows a small dip 
in the reststrahlenband of SiC, where the contribution of 
the evanescent modes is already relatively large. The sur- 
face roughness correction is in this case extremely small 
and originates in equal parts from the correction to the 
evanescent and propagating modes. 



0.04 



0.02 



AD E - 
D (scaled) 
D E pr (scaled) 






V 



0.1 0.5 



1.5 

n« -1 



2.5 



Figure 11: (Color online) Plot of D E and AD E as denned in 
Eq. (|10|) and Eq. (|34|1 for propagating and evanescent modes 
(dashed line) and for propagating modes only (dotted line) 
using <5 = 5 nm and a = 200 nm and distance z = 5 /im. 



For very small distances the LDOS is solely dominated 
by the evanescent contribution so that in this case the 
spectrum is quite different from that of the propagating 
par t 21 ' 24 . In Fig. [T2]we show a plot of D E for z = 500 nm, 
i.e., in the evanescent regime, leaving the other param- 
eters unchanged. As can be seen, the spectrum of D E 
has a resonance due to surface phonons. The curve of 
AD E shows a negative deviation of about 12 percent at 
the surface phonon resonance and a positive deviation of 
about 6 percent slightly below and above that resonance 
as can be seen in the inset. This behaviour is due to the 
scattering of surface phonons discussed in the previous 
section resulting in a broadened surface phonon disper- 
sion. 



A D 

D E (scaled) 




Figure 12: (Color online) Plot of D E and AD E as denned in 
Eq. (|34[) for z = 500 nm using the same roughness parameters 
as in Fig. fTT] 



D. Distance dependence of the LDOS 

Let us now focus on the distance dependence of the 
LDOS. For rough surfaces it can be expected that 
the roughness correction of the LDOS AD E defined in 
Eq. ([M]) converges for small distances, i.e., for z <C a, to 
the PA result. On the other hand, for large distances, 
i.e., for z 3> a, this correction should be inversely pro- 
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portional to the correlation length and can be described 
quantitatively utilizing the LDA of the proper self en- 
ergy in Eqs. ([22]) and ([23]). In Fig. [13] we show a plot 
of AD E for SiC at the frequency to = 10 14 s _1 using the 
roughness parameters 6 = 5 nm and a = 100 nm, 200 ran 
and 500 nm. It can be seen that the LDOS converges 
to the two limits for small and large distances giving a 
surface roughness correction on the order of some per- 
cent only for distances smaller than 100 nm. The arrow 
indicates that the surface roughness correction decreases 
when increasing the correlation length a. 




icr 8 icr 7 irr 6 io~ 5 

z (meters) 

Figure 13: (Color online) Plot of AD E over the distance for 
SiC using the frequency u = 10 14 s _1 for a rough surface with 
S = 5 nm for correlation length a = 100 nm, 200 nm, and 
500 nm. The thin solid line represents the PA result, whereas 
the thin dash-dotted line is the large distance approximation 
(LDA). 



PA 

a = 100nm 




icr 8 icr 7 io~ 6 icr 5 

z (meters) 

Figure 14: (Color online) As Fig. [Lj but for u = 1.787 • 
10 14 s" 1 . The red part of the curves indicates positive values 
and the blue one negative values. 

for u = 1.8- lO^s" 1 in Fig. [151 
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In Fig. [TJ] we show a similar plot when choosing the 
surface phonon frequency ui = uj s = 1.787 • 10 14 s _1 and 
plotting the modulus of AD E . As before, it can be ob- 
served that AD E converges to the PA and the LDA for 
z <C a or z 3> a, resp. Apart from that for distances 
z w a the surface roughness correction to the LDOS D E 
becomes negative, a feature not seen for tu = 10 14 s _1 . 
This behaviour can be understood from the above given 
discussion: the dispersion relation of the surface phonons 
is broadened due to roughness induced scattering yield- 
ing a negative AIm(r p ) for frequencies around the sur- 
face phonon frequency and for wavevectors around a -1 , 
because in that region the roughness induced scattering 
of surface phonons is strong. Since for a given distance 
z the main contributions to -D E stem from wavevectors 
k w z^ 1 , this roughness induced broadening becomes im- 
portant for distances z ~ a giving a smaller LDOS than 
for a flat surface so that AD E is negative. As can be seen 
in Fig. [14] the correction varies from about +10 percent 
to about —50 percent in the distance regime z > 100 nm 
for a = 100 nm. 

Now, as can for example be seen for AIm(r p ) in 
Fig. [9] (b) and for the spectrum of AD E in Fig. [12] for 
frequencies slightly below or above the surface resonance 
the LDOS increases due to the roughness induced scat- 
tering of surface phonons. Hence, relatively large positive 
surface roughness correction can be expected for such fre- 
quencies. To illustrate this statement we also plot AD E 



Figure 15: (Color online) As Fig. [13] but with u = 1.8 ■ 
10 14 s-\ 

It is apparent from the above given examples that the 
distance dependence of the surface roughness correction 
to the LDOS is not only sensitive to the surface roughness 
parameters themself, but also to the chosen frequency. 
Additionally, the surface roughness correction AD E is in 
general nonmonotonous and especially large for frequen- 
cies around the surface phonon resonance. 

VII. SUMMARY 

In this work, we studied the near-field radiative heat 
transfer between a nanoparticle and a rough surface uti- 
lizing direct perturbation theory up to the lowest nonva- 
nishing order in the surface profile. Employing the ma- 
terial properties of SiC we have shown that the distance 
dependence of the roughness correction to the heat flux 
is nonmonotonous and can be qualitatively understood 
from the roughness correction to the LDOS. 

We have derived an approximation in the small dis- 
tance regime, i.e., for distances d much smaller than the 
correlation length a, and have shown that it exactly co- 
incides with the results of the proximity approximation. 
Hence, the roughness correction is well described by the 
PA for d <C a. Therefrom, one can conclude that the 
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PA might also be helpful in other geometries as used in 
recent experimental setups^— to estimate the impact of 
surface roughness to the near-field radiative heat trans- 
fer at small distances. Since, the numerical results give 
larger values for the heat transfer than predicted by the 
PA for most distances d , one can expect to get an esti- 
mate of the lower limit of the surface roughness correction 
when using the PA. 

In the large distance limit d 3> a we have derived a 
simple approximation for the corresponding expressions 
of the LDOS and the heat flux, which simplifies the nu- 
merical calculations in this limit tremendously. We have 
shown that in this regime the rough surface can be re- 
placed by an equivalent surface layer of thickness S 2 /a for 
correlation lengths a smaller than the skindepth making 
contact with the results of Refsi 32 i 32 The corrections to 
the heat transfer are therefore relatively small in that 
regime when considering roughness with o>i. 

In the intermediate regime we could show that the 
LDOS and heat flux for a rough material can be smaller 
than that of a material with a fiat surface due to the 
roughness induced scattering of surface phonon polari- 
tons. Furthermore, we have shown that due to surface 
roughness the LDOS and therefore the heat flux has an 
s-polarized surface polariton contribution. The mecha- 
nisms behind these two unexpected results has been dis- 
cussed. 

Finally, we want to emphasize that the results for the 
LDOS presented in that work have a much larger range 
of applicability, since the LDOS, for instance, also deter- 
mines the lifetime of atoms and molecules near a surface. 



with 



G E(2 V) 



k 2 



1 - e 
~2tT 



dV 



x T(fc+ k+)K(k+, K')T(k-\ k-')f( K ', z') 
x 2( 7t ' - 7 ;)T(fc+, fc+)L(fc+, fc t -', kf) 



xT(K",k-")f(K",z') 



where 



27r 



S {n) {K-K') = Jd 2 xe Y - 



wherc i,j = (s,p) and 



a s (fc ± ) = z x k and a p (fe ) = a s 



Ik±| 



(A2) 



(A3) 

<"-"'> (S(x)) fl (A4) 



for n = 1, 2. Here all matrices A are defined by 



Aij (fcf , kf ) = Aij sa {kf ) <g) a,- (fcf ) , (A5) 



(A6) 
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are the normalized and orthogonal polarization vectors; 
z is the unit vector in z-direction, k = k/k and k^ = 

The matrix T is given by the diagonal matrix 



T 



' 2 7r 

7r+7t 







<57r+7t , 



t s 
t n 



(A7) 



Appendix A: Perturbation result for Green's 
function 



defining the amplitude transmission coefficients t s and 
t p . The components of the matrices ]N(/c I f ,fc t r ) and 
L(fcj + , fc t ~', fc<r") arc given by 



With the procedure in Refj^ 3 - and the perturbation the- 
ory in Ref.— it is possible to determine the Green's dyadic 
for each order with respect to the surface profile. The 
zcroth- and first order expressions can be found in Refi^ 3 -. 
For the second-order Green's dyadic we find for z < z' 



G k 



d 2 K 
(27T) 



^ e i(«.x+ 7r2 ) G E( 2 )^ ( M ) 



A ss = k ■ k'(j t +7^), 



N m 



{k X K%(7t +7t): 



1't 

(k x k ) z 2 , i\ 
(7re + 7t7t)> 



(A8) 
(A9) 

(A10) 



AW = -T2 [-ws'(7t + 7t'e) + « ' «'7t'(7 r 2 e + 7t7t')] • 

(All) 
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and 



Now. inserting the above defined components of the Ma- 
trices IN and IL yields 



L, 



_L sp — t=- — (k X K )z(& • K ) 

v efc o 



7t 



(k x k') 



lilt 111*1 A // 
K K + 7r7t K ' K 

_L 1 



x k ) z (k ■ k ) 



K ' 2 + 7 ; 7 ; 



7t'7? 



K K + K ■ K 7r7t 

7^1 «' 2 + 



(A12) 



(A14) 



(nn'y/e - 7t7^K ■ k') 



(A15) 



M s = -i(M) 2 (e-l)7t 

+ (fco<5) 2 (e-l) 2 fc 2 



fjj$B{\*-*'\) (B6) 



/J s °( K ')(^-«) 2 + 4^^^ xA ') 2 



(A13) M p = i(M) 2 ^-^7t( 1 



k e + 1 
« e 



- 7t 2 £> s ° («')(« x £') 2 + ^(k' + k • k 7^7t) 



(B7) 



fcge 



x (kk'c — k ■ k! 7t7t) 



Finally, these expressions can be further simplified by 
utilising cylindrical coordinates and introducing the 
modified Besscl functions 



W = ^-l d^cos(n^)e- xcos ^\ (B8) 



2tt 



Appendix B: Definition of the proper self-energy 



we get for s- and p-polarized modes the relation 

M s/p = M s/p , + M s/p>1 (B9) 



Inserting the second-order Green's function in Eq. (|A1|) 
into Eq. (j3J) and averaging allows for finding the second- 
order surface roughness correction to the LDOS in 
Eq. ([TU|) - (IT2l when defining the proper self-energy M s / p 



s/p 



(A„«) 2 (-i)(e-l) 



87 2 (^° /p ) 2 

(Bl) 



with iV^p 6 given by 



N s a (K }K ) = (t s ) 2 N ss (kt,K) 
N p l ( K ,K) = (t p ) 2 N pp (k+,K)h p 

iV p fc ( K , k', k) = (<p) 2 L P p(fc+, fc t "', K)h p 



(B2) 
(B3) 
(B4) 
(B5) 



with 



M s , = -i(M 2 7t(e- 1 ): (BIO) 

M Stl = eV dxxe~ x {kleD^xa^ 1 ) 
Jo 

x [Io(nax) + I2{kclx)\ (Bll) 
+ D°(2sa- 1 )7 1 .(2a;a- 1 )7 t (2a;a- 1 ) 
x [Jo (/cox) — ^(Kax)} }, 



and 



M Pi0 = i(fc 5) 2 7t- 



(e-1)/ K 2 e + 1 



1 - 



fc 2 6 



(B12) 



M p ,i = V I dMe- l2 | -j-^Dl(2xa 



e^-— j i (nax) 
-l . 



+ 7t(K)(7 r (2xa 1 )e — 7t(2xa 1 ))(2xa 1 K)/i(«aa;) 
- ^(Jo(kox) + 7 2 (Kas))7r(2a;a _:L )7t(2a;a~ :L )7 2 (K) 



D°(2xa _1 )7t (tt)[lo(Kaa:) - J 2 («ax)] 



(B13) 
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Appendix C: Approximation of the proper 
self-energy for \koy/ea/2\ <C 1 

Employing \kf)yJta/2\ <C 1 for 7 r ,7t, and we find 



grated analytically yicldings* 



for M Sj i given in Eq. (|B11I) 

where 



^-(iS(H+/ 2 °H)— -(i„ 2 H^H) 

(Cl) 



V = (k S) 



2 (e-l)2 

e 4 



(C2) 



and 



C(0 = / dxx m e- x I n (£x) 
Jo 

For the p-polarized modes we get 



2 e a 



V { ekfi e + 1 2 



(e- l)7t(K)- 5 K/ 1 2 (Ka) 



2k 



j 2 M + - 7t 2 W (/ 2 M + /|(ko)) 



W t {n)(l^a)-l" 2 {Ka)) 



(C4) 



The integrals of modified Bcsscl functions can be inte- 



e 
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(C5) 
(C6) 
(C7) 



By means of these relations one can implement both lim- 
its na -C 1 and ko > 1 by approximating the modified 
Bessel function. 



(C3) Appendix D: Large distance limit (z 2> a and a 2> d s ) 

With a similar procedure as above, but in the opposite 
limit |fcov / ^ a /2| > 1 we find for z a the lowest order 
terms of the self energy 



Af B ,i « (k 6) 



e 2 efc 2 



(-i 7t 2 fc ^), ( D1 ) 



:(^-l) 2 : 



ifc . 



(D2) 



As is apparent, in that limit the self energy is indepen- 
dent from the correlation length a. 



1 K. Joulain, J.-P. Mulet, F. Marquier, R. Carminati, and 
J.- J. Greffet, Surf. Sci. Rep. 57, 59 (2005). 

2 A. I. Volokitin and B. N. J. Persson, Rev. Mod. Phys. 79, 
1291 (2007). 

' ! E. A. Vinogradov and I. A. Dorofeyev, Physics Uspekhi 
52, 425 (2009). 

4 D. Polder and M. van Hove, Phys. Rev. B 4, 3303 (1971). 

5 L. Hu, A. Narayanaswamy, X. Chen, and G. Chen, Appl. 
Phys. Lett. 92, 133106 (2008). 

6 E. Rousseau, A. Siria, G. Jourdan, S. Volz, F. Comin, 
J. Chevrier, and J.-J. Greffet, Nature Photonics 3, 514 
(2009). 

7 A. Narayanaswamy, S. Shen, and G. Chen, Phys. Rev. B 
78, 115303 (2008). 

8 S. Shen, A. Narayanaswamy, and G. Chen, Nano Lett. 9, 
2909 (2009). 

9 R. S. DiMatteo, P. Greiff, S. L. Finberg, K. A. Young- 
Waithe, H. K. Choy, M. M. Masaki, and C. G. Fonstad, 
Appl. Phys. Lett. 79, 1894 (2001). 

A. Narayanaswamy and G. Chen, Appl. Phys. Lett. 82, 
3544 (2003). 

1 M. Laroche, R. Carminati, and J.-J. Greffet, J. Appl. Phys. 
100, 063704 (2006). 

2 M. Francoeur, M. P. Mengiig, and R. Vaillon, Appl. Phys. 
Lett. 93, 043109 (2008). 

3 U. F. Wischnath, J. Welker, M. Munzel, and A. Kittel, 
Rev. Sci. Instrum. 79, 073708 (2008). 



14 A. Kittel, U. F. Wischnath, J. Welker, O. Huth, F. Riiting, 
and S.-A. Biehs, Appl. Phys. Lett. 93, 193109 (2008). 

15 M. Nieto-Vesperinas, Scattering and Diffraction in Physi- 
cal Optics, (World Scientific, Singapore, 2006). 

16 A. A. Maradudin, Light Scattering and Nanoscale Surface 
Roughness, (Springer Science, New York, 2007). 

17 J.-P. Mulet, K. Joulain, R. Carminati, and J.-J. Greffet, 
Appl. Phys. Lett. 78, 2931 (2001). 

18 F. Pincemin, A. Sentenac, and J.-J. Greffet, J. Opt. Soc. 
Am. A 11, 1117 (1994). 

19 L. D. Landau and E. M. Lifshitz, Statistical Physics Part 
2, (Butterworth, Oxford, 2002). 

20 P.-O. Chapuis, M. Laroche, S. Volz, and J.-J. Greffet, 
Phys. Rev. B 77, 125402 (2008). 

21 K. Joulain, R. Carminati, J.-P. Mulet, and J.-J. Greffet, 
Phys. Rev. B 68, 245405 (2003). 

22 J.-J. Greffet, Phys. Rev. B 37, 6436 (1988). 

23 C. Henkel and V. Sandoghdar, Opt. Commun. 158, 250 
(1998). 

24 A. V. Shchegrov, K. Joulain, R. Carminati, J.-J. Greffet, 
Phys. Rev. Lett. 85, 1548 (2000). 

25 I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Se- 
ries and Products, (Academic Press, San Diego, 2007). 

26 C. Genet, A. Lambrecht, P. Maia Neto, and S. Reynaud, 
Europhys. Lett. 62, 484 (2003). 

27 P. Maia Neto, A. Lambrecht, and S. Reynaud, Europhys. 
Lett. 69, 924 (2005). 



13 



28 


B. N. J. Persson, B. Lorenz, A. I. Volokitin, Eur. Phys. J. 


35 


A. A. Maradudin and D. L. Mills, Phys. Rev. B 11, 1392 


20 


31, 3 (2010). 


36 


(1975). 


C. Henkel, Eur. Phys. J. D 35, 59 (2005). 


A. Marvin, F. Toigo, and V. Celli, Phys. Rev. B 11, 2777 


30 


G. Nogues, C. Roux, T. Nirrengarten, A. Lupascu, A. Em- 




(1975). 




mert, M. Brune, J.-M. Raimond, S. Haroche, B. Placais 


:S7 


S. J. Elston, G. P. Bryan-Brown, and J. R. Sambles, Phys. 




and J.-J. Greffet, EPL 82, 13002 (2009). 




Rev. B 44, 6393 (1991). 


31 


D. L. Brundrett, E. N. Glytsis, and T. K. Gaylord, Appl. 


38 


F. Marquier, C. Arnold, M. Laroche, J.-J. Greffet, and Y. 


:S2 


Opt. 33, 2695 (1994). 


39 


Chen, Opt. Exp. 16, 5305 (2008). 


T. S. Rahman and A. A. Maradudin, Phys. Rev. B 21, 504 


G. A. Farias and A. A. Maradudin, Phys. Rev. B 28, 5675 


33 


(1980). 


40 


(1983). 


T. S. Rahman and A. A. Maradudin, Phys. Rev. B 21, 


L. A. Moraga and R. Labbe, Phys. Rev. B 41, 10221 


34 


2137 (1980). 




(1990). 


G. S. Agarwal, Phys. Rev. B 15, 2371 (1977). 







